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We investigate cosmological constraints on an energy density contribution of elastic dark matter 
self-interactions characterized by the mass of the exchange particle msi and coupling constant asi- 
Because of the expansion behaviour in a Robertson- Walker metric we investigate self-interacting 
dark matter that is warm in the case of thermal relics. The scaling behaviour of dark matter self- 
interaction energy density (qsi oc a -6 ) shows that it can be the dominant contribution (only) in the 
very early universe. Thus its impact on primordial nucleosynthesis is used to restrict the interaction 
strength msi/ v / asi, which we find to be at least as strong as the strong interaction. Furthermore 
we explore dark matter decoupling in a self-interaction dominated universe, which is done for the 
self-interacting warm dark matter as well as for collisionless cold dark matter in a two component 
scenario. We find that strong dark matter self-interactions do not contradict super-weak inelastic 
interactions between self-interacting dark matter and baryonic matter (erf 10111 ^ ""weak) and that the 
natural scale of collisionless cold dark matter decoupling exceeds the weak scale (o"a M > °"wcak) 
and depends linearly on the particle mass. Finally structure formation analysis reveals a linear 
growing solution during self- interaction domination (S oc a); however, only non-cosmological scales 
are enhanced. 
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I. INTRODUCTION 

In the past decades high-precision observations al- 
lowed the development of a standard model of cosmology: 
ACDM. Its main statements are that we are living in a 
flat universe (Q? ot = 1.0052 ± 0.0064), dominated by the 
'dark' components: dark energy (f^ E = 0.721 ± 0.015) 
and non-baryonic dark matter (^q M = 0.233 ± 0.013) 

Q. 

The necessity of a dark energy component comes from 
the acceleration of the universe expansion, inferred from 
a high rcdshift Hubble diagram of type la supernovae as 
standard candles [|[ and radio galaxies as standard yard- 
sticks [|[ . A recent and impressive proof for the existence 
of dark matter (DM) can be deduced from observations 
of colliding galaxy clusters. Optical and near infrared 
observations of the galaxies, X-ray emission of the up- 
heated intergalactic plasma and gravitational lensing of 
the mass distribution show the necessity of a non- visible 
DM component, which dominates the mass budget [HQ. 
An overview of DM physics andparticle candidates can 
e.g. be found in recent reviews p-Q- 

Numerical structure formation simulations in the 
ACDM framework (from one of the first and most pop- 
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ular [l(| to the most recent [IT|, Ell ) show an impressive 
agreement with observations and have therefore become 
a cornerstone of modern cosmology. Nevertheless they 
also reveal two shortcomings that are worth being taken 
seriously. First, simulations predict scale-independently 
a large number of substructures in collisionless cold dark 
matter (CDM) halos, which exceed on galactic scales 
clearly the number of yet observed Milky Way satel- 
lites [III, [IH . One explanation is that reionization could 
prevent formation of visible baryonic structures in the 
smallest DM halos (e.g. [HI)- Second, simulations show 
cusps in the center of collisionless CDM halo density pro- 
files. But observations of dwarf spheroidal galaxies - 
which have a huge mass-to-light ratio and hence are ob- 
jects suited to study DM properties without perturbing 
baryonic effects, rotation curves of high spatial resolu- 
tion and large extension of low luminosity spiral galax- 
ies, and the universal rotation curve for spiral galaxies 
indicate a constant DM halo core density (e.g. fl^ - fl8l ]). 
An overview about processes that might lead from intrin- 
sic cuspy CDM distributions to the observed cored ones 
gives e.g. Ref. [r|. 

An idea to avoid both mismatches of the CDM sce- 
nario is to introduce strong elastic DM sclf-intcractions 
[20j | . A recent overview concerning collisional DM is 
given in Rcf. Q. Here, we want to concentrate on 
the most important facts that are also relevant for this 
work. The original self-interaction strength proposed 
by Ref. [H is cr S i/m D M = 0.45 - 450cm 2 /g (self- 
interaction cross-section over DM particle mass). But 
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in Rcf. [2l| it was demonstrated that cross-sections gen- 
erating reasonable dwarf galaxy cores predict too large 
galaxy cluster cores. Ref. [221 showed that indepen- 
dent of the dependence on the halo velocity dispersion 
self-interacting cross-sections cannot solve the satellite 
problem accurately. The most reliable constraint on the 
self-interaction strength comes again from observations 
of colliding galaxy clusters [23], [24[ ■ The nonexistence of 
an offset between the galaxy distribution and the gravi- 
tational lens mass peak, and the subcluster mass-to-light 
ratio allow to constrain osi/todm < 0.7cm 2 /g. This re- 
sult nearly completely rules out the formerly proposed 
self-interaction strength. The strongest limit on the col- 
lisional character of DM (asi/mnu < 0.02cm 2 /g) can 
be inferred from the observed ellipticity of DM halos and 
the property of DM collisions to make halos spherical [25| , 
but one has to take into account its model dependence 

[HSU. 

Another approach that avoids the satellite and cuspy 
halo problems is to provide the DM particles with a fi- 
nite thermal streaming velocity, to achieve a cut-off in the 
power spectrum and smearing of the innermost, highest 
density halo regions (2(|[27j (see also [HI). This means to 
introduce warm dark matter (WDM) particles. A lower 
bound on the DM particle mass can be determined from 
the Lyman-a forest (mwDM ^ 4keV, (2^|) and gravita- 
tional lensing (mwDM 2.2 kcV, [2Jj]) of high redshift 
quasars (given limits are for thermal relics). Ref. (30j 
showed that these boundaries can be lowered consider- 
ably in a set-up of mixed cold and warm DM, which we 
also consider in this work. So we follow Ref. |31[ in using 
1 — 10 keV as a typical mass range for WDM particles in 
the following. 

Other solutions proposed are stronger CDM annihi- 
lations ([HI, but see also [33j]) or a coupling between 
quintessence dark energy and CDM (e.g. |34j). 

Present attempts to enlarge the phenomenology of DM 
physics are e.g. strong DM baryonic matter interactions 
[35| or a dark radiation (electromagnetism) between DM 
particles [Hj]. 

We introduce in this work an energy density contribu- 
tion of elastic dark matter self-interactions. Despite the 
fact that self-interacting dark matter (SIDM) may not 
solve the shortcomings of the collisionless approach, the 
motivation for this work is to explore new, interesting 
cosmological consequences of an additional energy den- 
sity contribution of DM self- interactions within the above 
mentioned constraints. Interestingly enough, an interac- 
tion strength of usi/m ~ 1 cm /g still corresponds to 
strong interactions (<7 s t ron g ~ 10 fm 2 ) between nucleon- 
like particles (m ~ 1 GcV) . 

In Sec. |TT] we introduce our idea of a self-interaction en- 
ergy density contribution psi- Energy density scaling ac- 
cording to the Friedmann equations and its equation of 
state (psi = £?si) imply that the self- interaction contri- 
bution shows the steepest decrease with the scale factor 
(qsi oc a -6 ) and thus can (only) have a direct impact 
on the very early universe. Its proportionality to the 



SIDM particle density (gsi oc ng 1DM ) leads us to consider 
warm self-interacting dark matter (WSIDM) in the case 
of thermal relics to have the correct scaling behaviour 
("sidm oc a~ 3 ). But this does not rule out a second col- 
lisionless CDM component. 

After finally defining our parameter set, we use in Sec. 
IIIII today's DM energy density to constrain the pa- 
rameters characterising the SIDM particle properties and 
primordial nucleosynthesis limits on an additional en- 
ergy density contribution to constrain the self-interaction 
strength msi/y/otsi- We find that it depends inversely 
on the SIDM particle mass (msi/y/ a si <x 1/ w wdm) 
but can be at least as strong as for strong interactions 
(msiA/asi ~ 100 MeV). 

In Sec. IIVI we analyse the consequences on DM decou- 
pling in a universe dominated by the self-interaction en- 
ergy density contribution. The annihilation cross-section 
of WSIDM cr)^ DM is inverse proportional to the elastic 
self-interaction strength (cr)^ 01 ^ oc y^si/msi) and rather 
low (a^ DM <C Cweak) while the natural scale for the anni- 
hilation cross-section of a collisionless CDM component 
<7^ DM exceeds the weak scale and depends beside the 
self-interaction strength also on the particle mass tticdm ■ 
This casts new light on the 'WIMP miracle' and coincides 
with Fcrmi-LAT and PAMELA data (e.g. [Hill). We 
use the unitary bound and neutrino induced constraints 
on the DM annihilation cross-section to again limit the 
self-interaction strength. 

Another consequence of an early self-interaction domi- 
nated epoch may concern structure formation. We show 
in Sec. [V] that a rclativistic analysis of linear perturba- 
tion theory reveals a linear growing solution 8 oc a of 
self-interaction dominated SIWDM and also of collision- 
less CDM in a mixed model during self-interaction dom- 
ination. However, only non-cosmological scales (M < 
1O _3 M0) can be enhanced and a small observable effect 
could only be present with fine-tuned parameters. 
Finally we summarize our results in Sec. IVII 

II. SELF-INTERACTION ENERGY DENSITY 

In the following, we analyse constraints and conse- 
quences of an energy density contribution from dark mat- 
ter self-interactions qsi- 

We describe two particle interactions between scalar 
bosons or fermions by the exchange of vector mesons. 
For a scalar field 4> (fermionic field tp) and a vector field 
Vfj, the Lagrangian reads 

£ b = D;m-mb«-l^^ + | m v^- (la) 
Ct = $ (if -m)il>- jV^V^ + ImlV^V 1 * , (lb) 

with V^v = d^Vy — dpVft. The boson (fermion) field is 
coupled to the vector field by a minimal coupling scheme 

£V = <9 M + ig-v^Vn , (2) 

where <7 V 0(v>) is the <j>(i())-V coupling strength. We treat 
the vector field as a classical field. In homogeneous and 
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isotropic matter the spatial components of the vector 
field vanish and the equation of motion for the scalar 
(fermion) field reads: 



[V^W + m 2 h ] <j>{x) = 
[if - mt] ip(x) = 



(3a) 
(3b) 



In the mean-field approximation, after expanding into 
plane waves, we obtain the dispersion relation: 



k 2 



(4) 



Note that the vector interaction between the scalar 
(fcrmionic) particles is repulsive which ensures the over- 
all stability of selfinteracting boson matter and avoids an 
enhancement of the annihilation cross-section due to the 
formation of bound states. 
The number density of bosons (fermions) 

n h = J Q = i(j>*do<j> — i(d <l>*)<p-2gy < j,Vo<l>*<t>, (5a) 
m = J = V>7oV> , (5b) 

is just the source term for the vector field that is deter- 
mined from the equation: 

mlVo = g vcf> {i(j)*dQ(j)-i{d (j)*)(j)-2g v< pVQ(j)*(j)) 

= 9v<f> n b (6a) 
m v Vo = g V 4,^joi> = g v ^, n { (6b) 

The total energy density of the boson (fermion) matter 
can be determined from the energy-momentum tensor 



Qb - Qb +o m v l/ o — Qb + 



1 



5 V 
2m?, 



Qi = ^>7° (*<9 - ffv^Vb) i/j + -m*V$ 

Q 2 



(7a) 



(7b) 



where the equation of motion for the vector field (JSJ) has 
been used. The energy density of free boson matter obeys 
for the lowest energy mode k = 0, g h Tee = 2m 2 ) (j)*(f) = 
m h n h . 

The pressure is given by: 



Pb = Pb^ + 2 m t V ° = Pb^ + 7^2 "b V™) 

Pi = -^P [7° (ido - g V ijVo) - TO f ] ip + -m^ 2 

= \i> [7° Wo - g^Vo) -m^ + ^nj (8b) 

The pressure of free boson (fermion) matter fulfills p lree = 
g free /3 while the particles are relativistic, and for the 
lowest energy mode k = the total pressure is just given 
by the vector field contribution. 

The form of the interaction is equal to the one used for 
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investigating implications of interactions of fermions and 
bosons on compact stars in Refs. [39ll40j. 

For simplicity we denote the particle masses tosidm = 

^b(f)) ^SI — m vi 

and define the coupling constant asi = 
5v0(v>)/^' so that the energy density contribution from 
dark matter self-interactions reads 



QSI 



"SI 



SI 



'SIDM 



PSI 



(9) 



with msi/y/asi as the energy scale of the self-interaction. 
This scale can also be interpreted as the vacuum expecta- 
tion value of the Higgs field of the interaction. For weak 
interactions the interaction strength is given by Fermi's 
constant 1 fy/Gp or the vacuum expectation of the Higgs 
field generating the masses of the W and Z bosons, so 
that mweak/V^wcak ~ 300 GeV. Correspondingly, the 
strength of low energy strong interactions, quantum chro- 
modynamics (QCD), is controlled by the pion decay con- 
stant as 1/ f 2 in chiral perturbation theory, with f„ being 
the vacuum expectation value of the sigma field. Hence, 
for strong interactions the typical interaction energy scale 

IS ^strong /^strong ~ 100 Me V [H]. 

Note that according to the underlying propagator Eq. 
© is valid only when msi > 5Tsidm, so when the 
vector meson is non-relativistic. Otherwise the self- 
interaction energy density contribution scales like radi- 
ation (g SI oc n§ roM /T,f IDM oc Ig roM ). Once msi/y/asi is 
fixed, this can also be used as a boundary condition on 
the coupling constant for given Tsidm- A corresponding 
discussion is given in Appendix [XJ Please also note that 
the coupling constant asi defined here differs from the 
common definition by a factor 27r, e.g. a s = g 2 /(4ir). 
The equation of state ([9]) of the self-interaction psi = Qsi 
represents the stiffest possible equation of state consis- 
tent with causality. 

The expansion behaviour of the universe in a 
Robertson- Walker metric is described by the Fricdmann 
equations 



df? = g + p 
da a 



(10) 



where a is the scale factor, implying that the self- 
interaction energy density contribution shows the steep- 
est decrease (Fig.[T|): 



Qsi oc a 



-6 



(11) 



So, the universe could be in a self-interaction dominated 
epoch prior to radiation domination in the very early uni- 
verse. Under certain assumptions it might be possible in 
the future to constrain the dominating equation of state 
in the early universe via gravitational waves [4l| . 

The scaling behaviour gsj oc a~ 6 and the proportion- 
ality between self-interaction energy density and SIDM 
particle density gsi oc Hg IDM imply that: 

"Sidm oc a~ 3 (12) 

This condition is naturally fulfilled after DM decoupling 
and in the case of relativistic particles also before that, as 
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FIG. 1: Double- logarithmic plot of the evolution of different 
energy density contributions with the scale factor a. bcmb is 
the scale factor at photon decoupling, a£q dmat at radiation- 
matter equality and af. D . at the freeze-out of the neutron to 
proton number ratio, gsi is fixed so that self-interaction- 
radiation equality is at freeze-out of the neutron to proton 
number ratio. 




o (a ) 



FIG. 2: Evolution of the SIDM energy density £>sidm = 
Qsi + jwdm with the expansion of the universe for differ- 
ent self-interaction strengths. WDM particle parameters are: 
wwdm = IkeV, jwdm = 2, /avdm/2\vdm = 0, j^* c = 
1107, -Fwdm = 1- a t.o. is the scale factor at the freeze-out 
of the neutron to proton number ratio and ocmb at photon 
decoupling. 



"wdm oc T 3 . But it is incompatible with the exponen- 
tial suppression of non-rclativistic particles before they 
decouple (ticdm oc exp (— tticbm/T)). This is why we 
discuss in the following warm self-interacting dark matter 
and consider cold dark matter as a second, collisionless 
component. 

We describe the DM as a free Boltzmann gas and ad- 
ditionally take into account a non-zero DM chemical po- 
tential. Fermionic DM with a non-vanishing chemical 
potential according to Fermi-Dirac statistics is analyzed 
in detail in Ref. |42|]. Altogether our model contains six 
parameters. Besides the particle mass towdm, the degen- 
eracy factor <?wdm and the relativistic chemical potential 
^wdm/Jwdm, the number of degrees of freedom of parti- 
cles in thermal equilibrium at WSIDM decoupling g^q 
(which fixes today's WDM temperature) characterizes 
the SIDM particle properties. Today's relative amount 
of WSIDM is given by F^ DM = Q WDM /Q£ M , and finally 
m Si/V a si determines the WDM self-interaction strength. 



The evolution of the total SIDM energy density 
£?sidm = Qsi + £>wdm is shown in Fig. [2] Even very 
strong DM self-interactions can only have a direct impact 
on the very early universe. The cosmic microwave back- 
ground radiation (CMB) is unaffected by the additional 
DM self-interaction energy density contribution consid- 
ered. Primordial nucleosynthesis is the major corner- 
stone to constrain the self-interaction strength. Besides 
the scaling £>sidm oc a~ 6 while the self-interaction con- 
tribution dominates over the particle contribution, one 
sees that psidm oc a~ 4 once the relativistic particle con- 
tribution dominates and psidm oc a -3 when the WDM 
particles have become non-rclativistic. 



III. CONSTRAINTS 

A. Today's Dark Matter density ^dm 

The parameters describing the DM particle properties 
are not all independent of each other, and today's relative 
DM energy density can be used to extract allowed 
combinations. At the present-day temperature the WDM 
particles are non-relativistic, so that their energy density 
is given by £> WDM = m WDM ™ W dm- According to particle 
number conservation after decoupling and entropy con- 
servation of the particles in thermal equilibrium, today's 
WDM particle density depends on the one at decoupling 
as: 



'WDM 



, 10 

'TT 



„Wdc ( 
»th oq 



1 Q ^Wdcc 
If® n WDM 
-'Wdcc 



(13) 



Since we describe the WDM as a free Boltzmann gas, its 
number density until decoupling is given by: 

/r^x ffWDM t3 / MWDM \ , 1/n 

"wdm (T) = T cxp(^ - J (14) 

Thus Eq. (|13p leads to the following constraint on the 
WSIDM particle parameters 



5WDM m WDM / /AYDM \ 37T 

,Wdcc tpO CX Pl 



H 2 O 
2 "0 "DM 



3th oq WDM 



'VTwdm7 8x3if Wpi T 3 



O h 2 
1.80 eVx » , (15) 
0.1143 v ; 



with T = (2.725 ± 0.002) K [£| and DM /ig = 0.1143 ± 
0.0034 [l|. Allowed ranges of the parameters describ- 
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FIG. 3: Allowed combinations of the parameters describing 
the WSIDM particle properties according to Eq. (|15|1 . so re- 
sulting in the correct today's DM energy density Ddm — 
0.233. 



ing the WSIDM particle properties are shown in Fig. [3] 
where the degrees of freedom in thermal equilibrium at 
WSIDM decoupling is plotted against the WDM parti- 
cle mass over the relative amount of WDM . For the 
largest possible number of degrees of freedom in thermal 
equilibrium at WSIDM decoupling in the particle physics 
standard model, namely = 106.75, the WDM par- 

ticle mass is restricted to #wdm wwdm < 193 eV. To 
achieve WDM masses of 1 — 10 kcV, additional degrees of 
freedom at WDM decoupling are required. All in all this 
mass range demands at least 10 3 — 10 4 degrees of free- 
dom in thermal equilibrium at WSIDM decoupling. This 
is common to all WDM models and not specific to the 
one discussed here. Additional degrees of freedom appear 
e.g. in supersymmetric theories or in theories with extra 
dimensions or in string gas cosmology models (see e.g. 

El SI). 



B. Primordial nucleosynthesis 



Together with Eq. (fT5j) this can be transformed into a 
lower bound on the WDM particle mass: 

77i wdm . 135 10.75 4 / 3 a H§ n° DM fl t teq c ~ 1/3 n q \ 
— m Pi ^3 TCTt U»J 



WDM 



7tt 3 3±? 



T 3 AN V 



^22.6eVx i^M^ gth °q (19) 
0.1143 AN V v ' 

Ref. [1^ obtains as limit for the additional energy density 
during BBN (2a uncertainty): 



AN V < 0.3 



(20) 



The resulting constraints on the WSIDM particle param- 
eters are shown together with those from today's DM 
energy density in Fig. |4] While the lower right corner 
is excluded by today's DM energy density, the lower left 
one is ruled out by the allowed energy density (which 
scales like radiation) during BBN. Smaller WDM parti- 
cle masses require a higher temperature and thus a lower 
number of degrees of freedom in thermal equilibrium at 
decoupling to have the correct energy density today. In 
the grey shaded region at the left the increased WDM 
temperature would result in a too large energy density 
during BBN. 
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An additional energy density contribution at big bang 
nucleosynthesis (BBN) is parametrized as additional neu- 
trino families A7V„ 



BBN _ BBN 
t?rad — f? 7 



+ (3 + AN„) ^ BN 



BBN 



Qe 



(16) 



assuming that it scales like radiation (f) r ad <x a~ 4 )- In our 
model this is true for the energy density contribution of 
the WDM particles pwdm but the self-interaction energy 



density drops faster (ggi 



oc a 



J ). Hence, we use limits 



on AN V to constrain p W DM only (g^M ^ (? BBN ), 
resulting in the following constraint on the WDM particle 
parameters: 



ffWDM CX P 

„Wdcc 4 / 3 
Wtheq 



< 



7^ 4 



360 x 10.75 4 / 3 



AN„ 



(17) 



FIG. 4: Allowed combinations of the WSIDM particle param- 
eters according to Eqs. (|15|) and (|19[) . 

To constrain the DM self-interaction strength via BBN 
we have to focus on the earliest stage of BBN, which is the 
freeze-out of the neutron to proton number ratio. This 
occurs at temperatures of Tf. Q . ~ 800 keV. 
A larger energy density results in a higher expansion 
rate of the universe (H oc g 1 ^ 2 ). This implies that the 
equality between expansion rate and reaction rate of re- 
actions that keep the neutrons and protons in equilibrium 
is achieved earlier, i.e. at a higher temperature, and hence 
the neutron concentration at freeze-out is enhanced. The 
time when the deuterium bottleneck opens is not affected 
by an additional energy density of DM self-interactions 
(as gsi oc a~ 6 ). Hence, the period between freeze-out of 
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the relative neutron concentration (X n = n n j (n n + n p )) 
and the beginning of nucleosynthesis reactions is slightly 
extended (At = tbBBN — if.o.)- During this time span 
free neutrons decay. Nevertheless, also the neutron con- 
centration when the deuterium bottleneck opens is in- 
creased (A^ BBN = X£°-exp(-Ai/T n )). Nearly all neu- 
trons available for the nucleosynthesis processes are in- 
corporated into 4 He. For this reason, the primordial 
4 He mass fraction (Yp ~ 2X^ BBN ) is the ideal probe 
to constrain the self-interaction energy density contribu- 
tion and thus the self-interaction strength. The detailed 
calculation is given in Appendix [Bl 

We assume for simplicity that the freeze-out of the neu- 
tron to proton number ratio occurs in a radiation dom- 
inated universe, so that the DM self-interaction energy 
density contribution does not exceed the radiation con- 
tribution: 



a tot 



f.o. 

ftad 



= 1 



1 SI 



^rad' 



< x 



SI 



< 1, 

(21) 

where Xgf- = Qsi'/Qrad- According to the defini- 
tion of the DM self-interaction energy density, Eq. (f2"Tj) 
translates into the following constraint on the DM self- 
interaction strength: 



30 x 10.75 Tf.p ffwDM 




/ MWDM \ 
V ?WDM / 



f o -i/2 3V30 x 10.75 



SI 



rp E720° PO 

J-i.o fl o"DM WDM 



x 3 



1.70 x 10 6 eV 2 



feff 

,2 



T 3 
J 



^m^o Tf. /879kcV -^wdm 



m WDM 

n 



(22) 



0.1143 (x S i/0.279) 1/2 ™wdm 



For the second equality we have used the relation between 
the parameters describing the WSIDM particle proper- 
ties according to Eq. (fi"5|). 

A robust upper limit on the primordial 4 He abun- 
dance inferred from observations is (2a uncertainty, from 
Ref. E2): 



Yp < 0.255 



(23) 



This implies a constraint on the DM self-interaction en- 
ergy density contribution at the freeze-out of the neutron 
to proton number ratio (see Appendix [Bjl as: 



si 



< 0.279 



(24) 



The resulting constraint on the DM self-interaction 
strength according to Eq. (|2"2"|) is shown in Fig. [5] The 
grey shaded region is ruled out by the allowed addi- 
tional energy density contribution of DM self-interactions 
at the freeze-out of the neutron to proton number ra- 
tio. The limit on the self-interaction strength scales 
inverse proportional to the SIDM particle mass. Even 
an additional energy density contribution of DM self- 
interactions of the strength of the strong interaction 
(^strong A/^strong ~ 100 McV) is consistent with the pri- 
mordial element abundances. 
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FIG. 5: Constraint on the DM self-interaction strength ac- 
cording to Eqs. (|22[) and (|24[) by the permitted energy density 
at freeze-out of the neutron to proton number ratio as kick-off 
of the primordial nucleosynthesis. 



If we associate with the vector meson exchange interac- 
tion a cross-section given by 



osi 



m si 



(25) 



where s = 4i?g IDM in the center of momentum frame, 
with -Esidm ~ ?sidm (~ TOsidm) as the relativistic (non- 
rclativistic) single-particle energy, Fig. [5] shows the 'late 
universe' constraints on the dark matter collisional cross- 
section as discussed in Sec. HI comparable with the con- 
straint on the DM self-interaction strength from halo 
structure of Ref. [48|. Note that the constraints on 
fSi/^SiDM are valid only if all DM is self- interacting 
(i'wDM = 1), whereas our constraint on the DM self- 
interaction strength via primordial nucleosynthesis has 
a trivial dependence on the relative amount of SIWDM 
(msi/y^l oc 4 DM ). 



IV. DARK MATTER DECOUPLING 

Chemical decoupling occurs when the expansion rate 
of the universe exceeds the dark matter annihilation rate 



Ta = "DM (CTAV) 



(26) 



where (cjav) is the thermally averaged product of the 
total DM annihilation cross-section and the relative ve- 
locity of the annihilating DM particles. The expansion 
rate of the universe is determined by the dominant en- 
ergy density contribution (H oc g 1 ^ 2 ). In the standard 
model DM decoupling takes place in a radiation domi- 
nated universe. In an epoch prior to radiation domina- 
tion when the WDM self-interaction energy density con- 
tribution dominates (Fig. [1]), DM decoupling can occur 
during this self-interaction dominated era. 
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FIG. 6: Constraint on the DM self-interaction strength de- 
pending on the SIDM particle mass for two different rela- 
tive amounts of SIWDM according to Eqs (f22j) and pijl. 
together with the proposed ranges 20, 24, 25] of the DM colli- 
sional cross-section, for the case that all DM is self-interacting 
(Fwdm = 1) according to Eq. (f2"5")l 



A. Self-interacting Warm Dark Matter 

In the case of thermal WDM relics the DM particles 
are relativistic at decoupling and their annihilation rate 
is therefore Tf DM = n WDM aY DM . In a universe that is 
dominated by the WDM self-interaction energy density 
contribution, also the expansion rate is proportional to 
the SIWDM particle density (qsi oc n^, DM , Eq. ©), so 
that for Ta = H the WSIDM annihilation cross-section is 
independent on the particle parameters but determined 
by the elastic self-interaction strength: 



WDM 
CT A 




7.45 x 10" 



msi 
100 MeV 



(27) 



'■ ^"wcak 



This dependence of the WSIDM annihilation cross- 
section on the elastic self-interaction strength is shown 
in Fig. [7J The annihilation cross-section is given here in 
units of the cross-section for weak interactions, which is 
defined as: 



^wcak 



(28) 



« 3.18 x 10~ 12 GcV~ 2 « 1.24 x 10~ 39 cm 2 

One notes that the WSIDM annihilation cross-section is 
rather small for reasonable (elastic) interaction strengths 
(ct a vdm < o- wcak ). Hence, WSIDM decoupling in a self- 
interaction dominated universe reproduces naturally and 
consistently the 'super weak' inelastic coupling between 
the WSIDM and baryonic matter, required in Sec. IIII Al 
to obtain typical WDM particle masses of 1 — 10 keV. 



FIG. 7: Constraint on the WSIDM annihilation cross-section 
depending on the elastic self-interaction strength, in the case 
of WSIDM decoupling in a self-interaction dominated uni- 
verse, according to Eq. ()27|) . 



Our analysis of SIDM decoupling in a self-interaction 
dominated universe complies with the qualitative state- 
ment of Ref . |4!| "that the elastic scattering cross section 
cannot be arbitrarily small given a nonvanishing inelastic 
cross section" . 



B. Collisionless Cold Dark Matter 

For a CDM species the thermally averaged product of 
the total DM annihilation cross-section and the relative 
velocity between the annihilating DM particles is given 
according to Maxwell-Boltzmann statistics by: 



l„ „,\ _CDM 



\ m CDM J 



1/2 



(29) 



Hence, the conditional equation for the decoupling of 
collisionless CDM in a universe dominated by the self- 
interaction energy density contribution of SIWDM reads: 



CDM 



/ ™CDM \ 
V Tcdcc / 



-1/2 



,Cdcc 
WDM 



3 4mpi n^DM mgi 



(30) 



The WDM particle density is that of a decoupled, rela- 
tivistic free Boltzmann gas. If one inserts for the CDM 
particle density the particle density of a non-rclativistic 
particle species, one arrives at: 



mcDM 

Tcde 



cxp 



?Cdec / 



V3 7T 3 / 2 m P i H$n a BM 5th 



Cdcc 




ffCDM"*WDM TOSI 



If one takes for the CDM particle density the decoupling 
particle density according to its today's relic density (Eq. 
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TABLE I: DM parameter sets shown in Figs. [8] to [10] Com- 
pared to the reference set 1, we increase the CDM degeneracy 
factor <7cdm in set 2, the WSIDM particle mass ttiwdm in 
set 3 and the relative amount of WSIDM Jwdm m se t 4 and 
decrease the CDM particle mass bicdm in sets 5 and 6. 



Set 



mcDM 



3CDM 



HlWDM 



'WDM 



lOTeV 
lOTeV 
lOTeV 
lOTeV 
100 GeV 
lOGeV 



IkeV 
IkeV 
lOkeV 
IkeV 
IkeV 
IkeV 



0.1 
0.1 
0.1 
0.9 
0.1 
0.1 



T3| for the CDM component), one gets: 



CDM 
^wcak 



2 Hq ^DM ( TO CDM 

3 Wpi 



T 3 



1 



V 2cdec 
WDM ™WDM m S i 



1/2 



F° 

WDM 



(32) 



While in a radiation dominated universe the CDM an- 
nihilation cross-section depends logarithmically on the 
CDM particle mass, we find that er^ DM is proportional to 
mcDM when CDM decoupling occurs in a self-interaction 
dominated universe! Furthermore, a larger elastic self- 
interaction strength implies a higher expansion rate, so 
that the annihilation cross-section has to be larger, in 
order that the CDM fits its particle density today. 
One can solve Eqs. pip and (|32p iteratively for a chosen 
parameter set of DM particle parameters. The correla- 
tion between the annihilation cross-section of collisionlcss 
CDM and the self-interaction strength of elastic WDM 
self-interactions is shown in Fig. [5] 
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FIG. 8: Collisionless CDM annihilation cross-section in de- 
pendence of the elastic WDM self-interaction strength accord- 
ing to Eqs. (|31|l and (|32p for DM particle parameter sets given 
in Table U 

The value of msi/^/asi is bounded below by constraints 
from primordial nucleosynthesis according to Eqs. 



and (|24p . The upper limit of msi/y/oisi comes from re- 
quiring a self-interaction dominated universe at CDM de- 
coupling (ggj dcc > £rad CC )- This condition implies rather- 
strong elastic WDM self-interactions, which result in 
large collisionless CDM annihilation cross-sections, ex- 
ceeding considerably the cross-section of the weak scale. 
This is in contrast to what is called 'WIMP miracle' in 
the standard model, the fact that a CDM particle with 
a mass around 100 GeV fits to today's DM relic density 
with an inelastic cross-section of the weak scale. The 
nonobserved decay of a Z boson into two DM particles 
rules out CDM masses of mcDM ^ toz/2 ~ 45.6 GeV. 
The linear dependence of the CDM annihilation cross- 
section on the CDM particle mass can also be recognized 
in Fig. M (sec also Fig. [JTJ). The ratio of CDM parti- 
cle mass and temperature at CDM decoupling becomes 
slightly larger in a self-interaction dominated universe 
compared to CDM decoupling in a radiation dominated 
universe, depending on the CDM particle mass. We find 
values of mcDM /Tcdec between 25 — 35 for the parameter 
sets given in Table HI 

If CDM decoupling occurs in a radiation dominated 
universe, the natural scale of the velocity weighted 
mean annihilation cross-section is (u^v) ~ 3 x 
10~ 26 cm 3 s- 1 / (1 - F°, DM ) [H, For the decoupling 
of collisionless CDM in a universe dominated by the DM 
self-interaction energy density contribution this becomes: 



{<7 A V) = 



1/2 



'PI 



^WDM 
1 _ 



n^CDM 



2.77 x 10~ 23 cm 



WDM ™WDM TOSI 



(33) 



to C dm/10TcV lMeV 
mwDM/lkeV msi/y/asi 1 



F° 

WDM 

- F° 

WDM 



Hence, the degeneracy in the CDM particle mass is re- 
moved and the natural scale of the CDM annihilation 
cross-section depends also on the CDM particle mass. 
All in all the natural scale of CDM decoupling can be 
increased by some orders of magnitude when CDM de- 
coupling occurs in a self-interaction dominated universe, 
depending on the elastic WDM self-interaction strength 
(see Figs. [TTJl and ITT]) . Interestingly enough, such boosted 
CDM annihilation cross-sections are able to explain the 
high energy cosmic-ray electron-plus-positron spectrum 
measured by Fermi-LAT and the excess in the PAMELA 
data on the positron fraction (e.g. [13, [Hj], see Fig. ITT]). 
A general upper limit on the DM annihilation cross- 
section is set by the unitarity bound. For s-wave domi- 
nated annihilation this is (49. |51|: 



-in 



CDM 



(34) 



Together with the natural scale of the CDM annihi- 
lation cross-section, the unitary bound sets an upper 
limit on the CDM particle mass of thermal relics. For 
CDM decoupling in a radiation dominated universe this is 

^0 \V2 



tocdm < lOOTeVh - F\ 



WDM; 



If collisionless CDM 
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decoupling occurs in a self-interacting dominated uni- 
verse, the unitarity bound leads with the corresponding 
natural scale (Eq. (|33| ) to the following limit of the ther- 
mal relic CDM particle mass: 



mcDM < 



37T 1 



1/6 



1/3 ( "1CDM 



PI 



- F° 

WDM 

F° 

WDM 



V ^Cdec 
1/3 

1/3 



1/6 



'WDM 



m S i 



1/3 



(35) 



2.86 x 10 12 eV x 



/"si, 

^ IkeV lMeV 



/ m CDM \ 
V Tcdoc / 



1/6 



1 



WDM 

F° 
WDM 



1/3 



(36) 



The resulting upper limit on the thermal relic CDM par- 
ticle mass for CDM decoupling in a self-interaction dom- 
inated universe is shown in Fig. [§1 The maximum ther- 
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FIG. 9: Upper limit on the thermal relic CDM particle mass 
from the unitarity bound depending on the elastic SIWDM 
self-interaction strength for collisionless CDM decoupling in 
a self-interaction dominated universe for DM parameter sets 
given in Table U according to Eq. (1351) . 



tend to lead to less cuspy halo profiles, we use the Milky 
Way Halo Average neutrino constraint of Ref. [52[ to 
compare with the predicted CDM annihilation cross- 
sections when collisionless CDM decoupling occurs in a 
self-interaction dominated universe. Fig. [10] shows the 
velocity weighted mean CDM annihilation cross-section 
for decoupling in a self-interaction dominated universe 
according to Eq. (|33[) together with the unitarity bound 
and neutrino bound for given CDM particle mass. One 
realizes that the combination of superstrong WDM self- 
interaction strengths together with very heavy CDM par- 
ticle masses are ruled out. With regard to the neutrino 
bound of the DM annihilation cross-section one should 
consider that it assumes that all the DM is collisionless 
CDM, which is not the case in our model with a much 
more weakly annihilating WDM component, so that it 
sets very strong limits on the self-interaction strengths. 
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FIG. 10: Thermally averaged product of the total collision- 
less CDM annihilation cross-section and the relative veloc- 
ity between the annihilating CDM particles in dependence of 
the elastic WDM self-interaction strength for the DM parti- 
cle parameter sets given in Table HI according to Eq. (|33[) . 
Also shown are as upper limits the unitarity bound (ub) and 
neutrino bound (vb) for the chosen CDM particle masses. 



mal relic CDM particle mass for CDM decoupling in a 
self-interaction dominated universe depends on the elas- 
tic WDM self-interaction strength and on the DM parti- 
cle parameters. The limit on tocdm is well in the TeV 
range. 

Another way to constrain the annihilation of DM is via 
the appearance of thereby produced particles. Neutri- 
nos proved to be the most useful final state, since they 
provide a stringent but conservative upper limit on the 
DM annihilation cross-section independent of branching 
ratios [33|, [HJ . Assuming s-wave dominated CDM anni- 
hilation processes we can directly transfer the DM anni- 
hilation cross-section limits from today's neutrino signal 
to the early universe CDM decoupling. Since our model 
contains three ingredients, namely, warm dark matter, 
non-vanishing elastic dark matter self-interactions, and 
larger cold dark matter annihilation cross-sections, which 



Finally Fig. Qj] displays the velocity weighted mean 
CDM annihilation cross-section for decoupling in a self- 
interaction dominated universe according to Eq. (|33|) in 
the common presentation as a function of the collision- 
less CDM particle mass, together with the Halo Average 
neutrino constraint of Ref. [52| . the unitary bound ac- 
cording to Eq. (|34|) . the 2a contours for fits to Fermi and 
PAMELA data assuming annihilation only to of 
Ref. and the best-fit lines to the PAMELA data for 
annihilations to e + e~ and W + W~ , respectively, of Ref. 




V. STRUCTURE FORMATION 

Now we want to study the impact of a self-interaction- 
dominated epoch in the early universe on the evolution 



10 




FIG. 11: Thermally averaged product of the total collision- 
less CDM annihilation cross-section and the relative velocity 
between the annihilating CDM particles in dependence of the 
CDM particle mass for the DM particle parameter sets given 
in Table HU according to Eq. (|33|l . Also shown are as up- 
per limits the Halo Average neutrino bound of Ref. [12] (^b) 
and the unitarity bound according to Eq. (|34[) (ub), as well 
as the 2a contours for fits to Fermi (Fermi fi) and PAMELA 
(PAM fi) data assuming annihilation only to of Ref. [38|] 

and the best-fit lines to the PAMELA data for annihilations 
to e+e" (PAMe) and W + W~ (PAMW) of Ref. [H]. 



TABLE II: DM parameter sets shown in Fig. 1111 Compared 
to the reference set A, we increase the CDM degeneracy factor 
Scdm in set B, the WSIDM particle mass ttiwdm in set C, the 
relative amount of WSIDM -Fwdm in set D and the elastic 
WDM self-interaction strength in set E and decrease it in set 



Set 



rresi/ v 



ffCDM 



BlWDM 







A 
B 
C 
D 
E 
F 



lMeV 
lMeV 
lMeV 
lMeV 
IkeV 
100 MeV 



IkeV 
IkeV 
lOkeV 
IkeV 
IkeV 
IkeV 



0.1 
0.1 
0.1 
0.9 
0.1 
0.1 



of ideal fluid cosmological perturbations. We use uniform 
expansion gauge (UEG) that is free of unphysical gauge 
modes and has the two gauge invariant variables 5 and 
ip 1 which can be identified with the density contrast and 
a quantity related to the fluid velocity in the subhorizon 
limit (/c p h ^> H), respectively [54l [55j. 

The system of general relativistic evolution equations 
in UEG reads 



5[ 



3 (wi - 4 k . 3 (1 + Wi ) 

Si + -fiWi a. 37) 

a an a 

3w» - 1 ; 2 k (1 + Wj)k 

-^^- C ^ 6i -^H— a (38) 



l + 3d 



(*) +1(1 + -) 



-s, 



(39) 



where S is the density contrast, w — p/g is the equation 
of state, c s is the isentropic speed of sound, H is the con- 
formal Hubble parameter, k is the comoving wavenumbcr 
of the mode (k = fc p h a), and a is the perturbation of the 
lapse. Primes denote derivatives with respect to the scale 
parameter a. Eqs. (|37| and (|38|) apply to each decoupled 
fluid component i individually and the general relativis- 
tic analogue of the Poisson equation (|3"9")l connects them. 
In Eq. (f39|) the averaged quantities w = J2iPi/ Si Qii 

S = YsiiQiSij/YjiQi and c s = Y,i{ctiQ%Si)IY,i{QiSi) en- 
ter. 

For a single fluid with w = 1 , corresponding according to 
Eq. ([9J to self-interaction domination, the analytic solu- 
tion in the subhorizon limit (k/H 3> 1) is given by 

5 S i(a) oc a ■ (Acos(a 2 - 3tt/4) + Bsm(a 2 - 3tt/4)) , 

(40) 

i.e. an oscillation with linearly growing amplitude as 
shown by Ref. (55| . This is in contrast to density fluctu- 
ations in standard CDM that can only grow logarithmi- 
cally during radiation domination in the early universe. 

For self-interacting warm dark matter there are two 
relevant damping scales, collisional self-damping (sd) due 
to particle scattering at early times and free streaming 
(fs) at late times when the WDM elastic self-interaction 
rate Tsi has dropped below the Hubble rate. These two 
can be estimated via the expressions 



sd 



r S i(t)« 2 W 

icollapaa V^^f) dt 



a(t) 



(41) 
(42) 



as given, for example, in Refs. [56|,l57j. i co n apso denotes 
the time of gravitational collapse, and the rms velocity 
of the WDM particles is ^wdm = c = 1. Eq. ([iljl is 
only valid as long as i 8 d "C k s , i.e. as long as the particles 
can be treated as interacting. Any increase in the den- 
sity contrast in WSIDM produced at early times will be 
washed out at later times either due to collisional self- 
damping or due to an inevitable phase of free streaming 
after self-decoupling (sdec). 

In a mixed model of SIWDM and collisionless CDM 
the picture can differ because the CDM component allows 
some increase in density fluctuations to be stored. First 
of all let us examine the solution to the Eqs. (|5T|) to 
(|3"9"|) for a subdominant CDM component in a SIWDM 



background (tocdm = C 



s CDM 



0,w 



Wsi 



l,cf 



1): 



CDM 



a 



k 

^ttV'cdm 
an 

"0CDM 

a 

6 

"(£) 2 + 



3 

— a 
a 

k 
<5si 



a 



(43) 
(44) 
(45) 



The transition from superhorizon to subhorizon behav- 
ior happens very quickly in a SIDM background since 
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k/H = kph/H = (a/a™) , where a' fc n is the scale param- 
eter at horizon entry. The terms proportional to a can 
be dropped in the subhorizon limit (k/H » 1) and the 
equations simplify further: 

^CDM = -^2^cdm (46) 
a k 

^CDM = V'CDM (47) 

a 

The solution for "0cdm is easily found to be V>cdm = C/a 
with C being a constant. This then automatically yields 
the solution to Eq. (|4T5)) 

<5cdm = a ■ (c/af) + D , (48) 

with D being a constant. Interestingly, this means that 
subhorizon collisionless CDM density fluctuations will 
also grow linearly during a SIWDM dominated phase in 
contrast to a radiation dominated phase. Thus there will 
be a region of enhanced fluctuations at low masses in the 
matter power spectrum between the comoving wavenum- 
ber that is equal to the Hubble scale H eq at SIDM- 
radiation equality and the wavenumber that corresponds 
to the collisional self-damping scale k°2 of SIDM at the 
same moment. A quick estimate then yields the follow- 
ing CDM transfer function T(k) = A cq (k)/A [n {k) i.e. the 
ratio of the amplitude of the CDM density fluctuation 
A with wavenumber k at SIDM-radiation equality (eq) 
normalized to the amplitude at horizon crossing (in): 

T ( fc ) = \/5' fc sd q >fc>^ Cq (49) 

This results from the fact that in a SIDM background 
modes become subhorizon more quickly since k/H oc a 2 
while each subhorizon CDM mode only grows as Jcdm oc 
a. So while each subhorizon mode has increased in am- 
plitude by one order of magnitude, two additional orders 
of magnitude in wavenumber have become subhorizon. 
Therefore the spectrum is less steep than one might ex- 
pect. 

In Refs. [H, limits on the abundance of planet 
size dark matter objects in the galactic halo via grav- 
itational lensing from the MACHO and EROS surveys 
are given. These surveys are sensitive to dark matter ob- 
jects down to ~ 10~'Mq. In Fig. H2lwe show the affected 
mass range of collisionless CDM density fluctuations as 
a function of the WDM elastic self-interaction strength 
TOsi/V a si for two different values of the SIWDM parti- 
cle mass rayvDM = IkeV, 100 keV at ^wdm/Twdm = 
and -FyyrjM = 0-1. Also shown is the sensitivity limit 
of the MACHO and EROS surveys. The shaded areas 
are limited to the left by the requirement a^ rad < Obbn 
(Eqs. (|22p and ([24)) ) which also limits the largest struc- 
tures that can be affected to ~ 1.4 x 1O~ 3 M . Using 
Eq. (|4"TJ)) this means in turn that fluctuations on scales 
of 10~ 7 Af© are only enhanced by at most a factor of 



- (10" 3 /10" 7 ) 1/6 w 5. Therefore it is unlikely that an 
observable overproduction of planet sized objects could 
result from the proposed scenario. For most of the pa- 
rameter space the effect is far from being observable with 
present small scale dark matter surveys. In the limit of 
very weak WDM self-interactions the self-damping length 
becomes very large and WSIDM can become free stream- 
ing before SIDM-radiation equality and any temporary 
small scale increase in the density contrast gets damped 
away. 

Note that elastic scattering processes between collision- 
less CDM and standard model particles can contribute an 
additional induced collisional damping scale until CDM 
thermal decoupling, discussed e.g. in Refs. [54l l56l [60| . 
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FIG. 12: Mass range of collisionless CDM density fluctuations 
affected by a SIWDM dominated epoch for two values of the 
SIWDM dark matter mass as a function of the self-interaction 
strength at ^wdm/7wdm = and F-^dm = 0.1. 



VI. CONCLUSIONS 

In this paper we have analysed constraints on an en- 
ergy density contribution of elastic dark matter self- 
interactions gsi, characterized by the mass of the ex- 
changed particle msi and the coupling constant cvsi- 
The scaling of energy densities implied that the self- 
interaction contribution decreases as gsi oc a -6 and thus 
can only have a direct impact on the very early universe. 
As the energy density scales with the number density 
squared due to interactions, self-interacting dark matter 
has to be warm in the case of thermal relics to give the 
correct scaling behaviour nsiDM oc a~ 3 . Note that this 
does not rule out a second collisionless cold dark matter 
component. 

We used today's dark matter energy density and the 
allowed radiation energy density during primordial nu- 
cleosynthesis to constrain the parameters characterising 
the warm self- interacting dark matter particle properties. 
The dependence of the primordial 4 He abundance on the 
dark matter self-interaction energy density contribution 
at neutron to proton number ratio freeze-out allowed to 
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constrain the self- interaction strength rnsi/y/oisi, which 
depends inversely on the self-interacting dark matter par- 
ticle mass (msi/y/cxsi oc 1/wwdm) but can be at least 
as strong as the strong interaction scale (msi/Vasi ~ 
100 Me V). Furthermore, our constraint on the dark mat- 
ter self-interaction strength has a trivial dependence on 
the relative amount of self-interacting warm dark matter 

(msiA/asT a F wdm)- 

We also analyzed dark matter decoupling in a uni- 
verse dominated by the self-interaction energy density 
contribution. The annihilation cross-section of warm self- 
interacting dark matter cr)^ DM is inverse proportional to 
the elastic self- interaction strength (er)^ DM oc v / ^si/ TO Si) 
and much smaller than a WO ak- The natural scale for the 
annihilation cross-section of a collisionless cold dark mat- 
ter component ct^ dm exceeds the weak scale (tr^ DM > 
Cwoak) and depends linearly on the particle mass mcDM 
(cr^ DM oc m CDM x y/a$i/mg\). This casts new light 
on the 'WIMP miracle' and coincides with the Fcrmi- 
LAT and PAMELA data. The unitary bound and neu- 
trino induced constraints on the dark matter annihila- 
tion cross-section allowed to disfavour the combination 
of superstrong elastic warm dark matter self-interactions 
(msiA/asi % 1 MeV) together with very heavy thermal 
relic cold dark matter particle masses (wcdm ~ 10 TeV). 
A relativistic analysis of linear perturbation theory re- 
veals a linear growing solution 8 oc a of self-interaction 
dominated warm dark matter and also of collision- 
less cold dark matter in a mixed model during self- 
interaction domination. However, only non-cosmological 
scales (M < 10 _3 M Q ) can be enhanced and a small ob- 
servable effect could only be present with fine-tuned pa- 
rameters. 



tosi > 5Tsidm- Once msi/^/asi is known, this can also 
be used as a boundary condition on the coupling constant 
for given T S idm: 

a SI > 25 ^ T| IDM (Al) 
m si 

This condition has to be fulfilled at primordial nucle- 
osynthesis and at collisionless CDM decoupling for the 
two component DM scenario. Fig. [T3] shows the corre- 
sponding constraints on the self-interaction coupling con- 
stant for the conservative assumptions T^p M = Tf . and 
t wdm = ™cdm/ (m CDM /T Cdcc ) with tocdm = 10 TeV. 
We sec that a self-interaction dominated universe at BBN 




m sl /<4| 2 (eV) 



FIG. 13: Constraints on the coupling constant in dependence 
of the DM self-interaction strength at primordial nucleosyn- 
thesis and collisionless CDM decoupling, according to Eq. 
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Appendix A: Self-interaction coupling constant 

The ansatz for the self-interaction energy density that 
we have used in this work (Eq. ^ is valid only when 



imposes rather loose conditions on the coupling con- 
stant (e.g. a s trong > 1-93 x 10~ 3 ), while collisionless 
CDM decoupling requires exotic coupling constants (e.g. 
a strong > 2.47 x 10 8 ). 



Appendix B: Primordial nucleosynthesis 

The following analytical calculations are based on Refs. 
[U [(32j . We assume that the freeze-out of the neutron 
to proton number ratio occurs in a radiation dominated 
universe, so that the DM self-interaction energy density 
contribution does not exceed the radiation contribution: 



f.o. 
ftot 



ft f a°d = (1 



' SI 



erad' 



o < x'gf < 1 



The modification of the total energy density implies a 
change in the temperature-time relation, and thus the 
conditional equation of the neutron to proton ratio freeze- 
out temperature becomes: 



9.50 



Tf.p. 
Am 



4.63 



Tf.p. 
Am 



0.677 



f.o.V2 
Scff 



Tt. . 
Am 



-f.o. 
^SI 



1/2 



(Bl) 
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gig' is the effective number of degrees of freedom con- 
tributing to the radiation energy density. According to 
the discussion in Sec. IIIIBI g f c ^- = 11.275. Adopting 
the largest possible energy density of radiation during 
BBN, we calculate the most conservative limits on the 
self-interaction energy density and hence self-interaction 
strength. The freeze-out temperature without any self- 
interaction energy density contribution - but with the 
maximum WDM particle contribution p^ DM = 0.3 gfc°- 
- is T f . .(4f = 0) ss 848 keV. Doubling the radia- 
tion energy density at neutron to proton number ratio 
freeze-out, the freeze-out temperature would increase to 
938 keV. The increase of the freeze-out temperature by a 
non-vanishing DM self-interaction energy density contri- 
bution at freeze-out results in an increase of the relative 
neutron concentration at freeze-out, which is given by: 




9.50ar 



(B2) 



For a vanishing energy density contribution of DM self- 
interactions the relative freeze-out neutron concentration 
is X£ ° (4f = 0) w 0.157 which would increase to 0.184 
if one doubles the radiation energy density. The num- 
ber of neutrons available for the primordial nucleosyn- 
thesis processes depends on the time spent between the 
freeze-out and the opening of the deuterium bottleneck. 
The moment of the neutron concentration freeze-out (us- 
ing the point in time corresponding to the above defined 
freeze-out temperature) in a radiation dominated uni- 
verse is: 



tf.o. 



45 

16 7T 3 



1/2 



J.O.-V2 



1 + X 



f.oA-V2 
SI ) 



T, 



r.o. 



(B3) 

The time corresponding to the temperature when the nu- 
cleosynthesis processes effectively set in is given by 



tbBBN 



45 

16 7T : 



1/2 



™P1 Scff 



bBBN- 1 / 2 rp-2 



bBBN 



(B4) 



since the energy density contribution of the self- 
interaction is vanishing by then (see the discussion in 
Sec. IIIIB"1) . Because T^bbn — 73.7 keV the correspond- 
ing effective number of rclativistic degrees of freedom is 



gbBBN ^ 3_5Q ; w here we include the WDM particles as 
being still rclativistic. For relatively heavy WDM parti- 
cles that decouple at a large number of degrees of free- 
dom in thermal equilibrium this does not necessarily be 
the case, but it is again the right choice for a conser- 
vative constraint of the self-interaction strength. With 
this input one can calculate the relative neutron concen- 
tration at the effective beginning of the nucleosynthesis 
( X bBBN = exp (-At/T n )) and finally the expected 
primordial abundance of 4 He (Y P ~ 2A^ BBN ). The de- 
pendence of the primordial 4 He abundance on an addi- 
tional energy density contribution of DM self- interactions 
at freeze-out of the number ratio of neutrons and pro- 
tons is shown in Fig. [TJ] Without this contribution it is 
Yp(xg°' = 0) « 0.241 and in the case of twice the radia- 
tion energy density this value increases to 0.281. 
The primordial 4 He abundance inferred from observa- 



0.28 



0.27 



0.26 



0.25 



0.24 



FIG. 14: Primordial 4 He mass abundance Y p depending on 
the additional energy density of DM self-interactions at neu- 
tron to proton number ratio freeze-out. Additionally the up- 
per limit on Y p inferred from observations is shown. 



tions is subject of systematic uncertainties (for a dis- 
cussion see Ref. [471) b ut a robust upper limit on Yp 
is Yp < 0.255 (2cr, |47j). This implies a constraint on 
the DM self-interaction energy density contribution at 
the freeze-out of the neutron to proton number ratio of 




r f .o. 



< 0.279. 



Very recently Refs. [63j, |64j determined the primordial 
4 He abundance Yp with a central value of 0.256, which 
underlines the possibility of new physics beyond standard 
BBN. 
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